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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

1 
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__.  ..  1- 

acre-feet 
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cubic  meters 

degrees  (angle) 
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acres 

j  nches 

25.4 

millimeters 

miles  (U.S.  statute) 

_ 

1.609347 

kilometers 

1  Introduction 


Background 

The  U.S.  Army  Engineer  Division,  Missouri  River  (MRD),  controls, 
maintains,  and  conserves  water  resources  on  the  mainstem  Missouri  River 
to  fulfill  project  purposes  authorized  from  1930  through  1940.  Since  au¬ 
thorization,  considerable  demographic,  social,  economic,  and  political 
changes  have  occurred  in  the  region.  In  1990,  MRD  began  re-evaluating 
the  Missouri  River  Master  Water  Control  Manual  to  identify  the  operating 
plan  that  best  meets  the  wide  variety  of  needs  in  the  basin.  Impact  assess¬ 
ment  methods  that  would  allow  MRD  to  identify  the  effects  of  different 
operating  plans  on  basin  resources  or  uses  had  to  be  developed.  These 
methods  would  facilitate  trade-off  analyses  to  help  MRD  identify  the  oper¬ 
ating  plans  that  would  provide  for  the  equitable  use  of  resources  for  au¬ 
thorized  purposes  (hydropower,  flood  control,  water  supply,  navigation, 
water  quality,  recreation,  and  fish  and  wildlife). 


Purpose 

The  authors  of  this  report  sought  to  develop  a  method  for  predicting  the 
impacts  of  system-operating  alternatives  on  fish  in  the  six  main  stem  reser¬ 
voirs  (Fort  Peck,  Sakakawea,  Oahe,  Sharpe,  Francis  Case,  and  Lewis  and 
Clark)  of  the  upper  Missouri  River. 


Objectives 

First,  our  aim  was  to  use  correlation  and  regression  to  quantify  the  ef¬ 
fects  of  seasonal  or  annual  variations  in  reservoir  hydrology  on  catches  of 
young-of-year  (YOY)  fish  in  summer.  Second,  we  hoped  to  develop  soft¬ 
ware  that  would  quickly  calculate  a  fish  reproduction  index  (RI)  for  every 
possible  year  in  the  period  of  record  (1898-1990)  for  any  operational 
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alternative.  We  wanted  the  MRD  to  be  able  to  evaluate  operational  alter¬ 
natives  by  comparing  a  long  chronology  of  predicted  indices  and  statistics. 
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2  Methods 


Hydrologic  Data 

We  developed  quadratic  relations  to  predict  reservoir  surface  area  and 
volume  from  the  surface  elevation  of  the  six  impoundments  (Table  1). 
These  relations  were  needed  because  end-of-month  data  on  area  and  vol¬ 
ume  were  not  as  readily  available  as  data  on  elevation,  inflow,  and  re¬ 
lease.  We  wanted  to  derive  independent  hydrologic  variables  from  area  or 
volume  instead  of  elevati  jn  so  that  their  dimensions  would  be  consistent 
with  area  or  volumetric  dimensions  associated  with  measures  cf  nutrient 
loading  and  reservoir  productivity.  Also,  fish  catch  per  unit  effort  in  gears 
like  seines  also  could  be  expressed  on  an  areal  basis. 

Hydrologic  data  consisting  of  end-of-month  elevations,  inflow,  release, 
and  subbasin  inflow  from  1967  to  1990  were  provided  for  every  reservoir 
by  the  MRD.  From  these  data,  we  derived  22  hydrologic  variables  describ¬ 
ing  annual  and  seasonal  hydrologic  characteristics  that  were  believed  to 
be  important  determinants  of  YOY  fish  catch  in  annual  samples  taken  by 
State  fishery  biologists  in  the  basin  (Table  2).  Plots  of  standard  deviations 
versus  means  of  all  independent  variables  indicated  which  variables  re¬ 
quired  transformation  to  stabilize  variances  and  normalize  distributions. 

We  checked  normality  using  the  UNIVARIATE  procedure  (SAS  Institute, 
Inc.  1988a).  Many  variables  were  transformed  by  taking  the  base  10  loga- 
iithm  of  values  plus  one.  Other  transformations  such  as  natural  logarithm, 
reciprocal,  and  square  root  were  tested  but  failed  to  provide  significantly 
better  normalization  than  the  base  10  logarithmic  transformation.  Change- 
in-area  variables  were  not  transformed  because  we  wanted  to  retain  both 
negative  and  positive  values. 

Hydrologic  variables  derived  for  each  reservoir  and  select  sample  statis¬ 
tics  are  presented  in  Tables  3-8.  We  explored  conditions  unique  to  individ¬ 
ual  reservoirs  to  account  for  differences  in  times  of  fish  spawning,  which 
varied  with  latitude,  by  using  different  sets  of  hydrologic  variables.  We 
did  not  use  area  variables  for  Lake  Sharpe  and  Lewis  and  Clark  Lake  (Ta¬ 
bles  6  and  8)  because  area  does  not  vary  significantly  among  months. 
Variables  for  these  riverlike  reservoirs  were  based  upon  inflow  and  flush¬ 
ing  rate  (total  release/mean  volume). 
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We  also  derived  independent  variables  similar  to  those  in  Table  2  from 
daily  hydrologic  data.  Our  goal  in  using  daily  data  was  to  determine 
whether  similar  models  would  be  derived  from  daily  and  monthly  data  and 
to  define  potential  limitations  of  using  end-of-month  data.  The  source  of 
hydrologic  data  for  evaluating  alternative  operations,  MRD’s  Long  Range 
Study  (LRS)  model,  provided  end-of-month  data  exclusively.  Two  vari¬ 
ables  unique  to  daily  data  included  maximum  change  in  area  and  the  coef¬ 
ficient  of  variation  (CV)  in  area  from  15  April  through  15  May.  We  hoped 
to  capture  negative  effects  of  short-term  drops  in  water  levels  that  might 
damage  fish  reproduction  by  disrupting  spawning  (June  1970,  Vogele 
1975,  Walburg  1976),  exposing  eggs  (Aass  1960;  Heman,  Campbell,  and 
Redmond  1969;  Priegel  1970;  Estes  1971),  or  concentrating  YOY  for  pred¬ 
ators  (Bennett  1962,  Jenkins  1970,  Beard  and  Snow  1970,  Aggus  1979). 


Weather  Data 


Temperature,  wind,  and  storm  frequency  are  believed  to  be  important 
factors  affecting  the  reproductive  success  of  many  reservoir  fishes  (Wal¬ 
burg  1972,  Clady  and  Hutchinson  1975,  Clady  1976,  Summerfelt  1975, 
Nelson  and  Walburg  1977,  Aggus  1979).  Wind  and  waves  can  increase  tur¬ 
bidity  and  sedimentation  along  shorelines,  and  sedimentation  adversely  af¬ 
fects  survival  of  eggs  and  YOY  fish  (Hassler  1970).  Weather  data 
recorded  hourly  from  1973  to  1990  at  seven  municipal  airports  a'ong  the 
upper  Missouri  River  were  obtained  from  the  National  Oceanic  and  Atmo¬ 
spheric  Administration,  U.S.  Department  of  Commerce.  Sites  included 
Glasgow,  MT,  for  Fork  Peck  Lake;  Williston  and  Bismarck,  ND,  for  Lake 
Sakakawea;  Bismarck,  ND,  and  Pierre,  SD,  for  Lake  Oahe;  Pierre,  SD,  for 
Lake  Sharpe;  and  Pierre,  SD,  and  Norfolk,  NE,  for  Francis  Case  and 
Lewis  and  Clark  Lakes.  We  used  a  heat-exchange  model  to  calculate  daily 
equilibrium  temperatures  from  estimates  of  percent  cloud  cover,  air  tem¬ 
perature  (wet-bulb,  dry-bulb,  and  dew-point),  wind  speed,  longitude,  lati¬ 
tude,  and  elevation.  Equilibrium  temperature  was  used  as  a  surrogate  for 
water  temperature,  which  was  rarely  available.  We  derived  three  indepen¬ 
dent  variables  including  equilibrium  temperature,  the  frequency  of  wind 
speeds  exceeding  the  75th  percentile  wind  speed,  z  id  the  number  of  storm 
hours  from  31  March  through  30  June  every  year,  coincidental  with  fish 
spawning  and  nursery  periods.  A  storm  hour  indicates  that  a  storm  was 
present  during  some  part  of  an  hour  monitored  by  a  weather  station.  Equi¬ 
librium  temperature  was  log  ^-transformed;  a  square  root  transform  was 
used  on  wind-speed  and  s tor" -hour  variables. 
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Covariates 


We  used  sport-fish  stocking  data  from  Montana,  North  Dakota,  and 
South  Dakota  to  derive  variables  for  use  as  covariates  in  regression  analy¬ 
sis.  Variables  included  base  10  logarithmic  transformations  of  one  plus 
the  number  of  fmgerling  or  fry  of  various  species  stocked  annually. 


Intercorrelations 


Because  of  the  nature  of  annual  and  seasons’  hydrologic  events,  most 
hydrologic  variables  exhibited  some  degree  of  intercorrelation.  When  in¬ 
dependent  variables  are  correlated,  regression  coefficients  are  not  unique, 
but  depend  on  other  intercorrelated  variables  in  the  model.  Nevertheless, 
correlations  between  independent  variables  usually  are  not  a  serious  prob¬ 
lem  if  the  goal  is  to  derive  models  for  inference  or  to  predict  new  observa¬ 
tions  (Neter  and  Wasserman  1974). 

Low  degrees  of  intercorrelation  were  accepted,  so  we  could  use  as 
many  variables  as  possible  for  regression.  We  used  intercorrelated  vari¬ 
ables  when  they  explained  less  than  SS  percent  of  the  variation  in  other  in¬ 
dependent  variables  (r  <  0.75;  r2  <  0.55,  where  r  =  correlation  coefficient, 
r2  =  coefficient  of  determination).  We  forced  regression  models  to  use 
only  one  of  several  more  highly  correlated  (r2  >  0.55)  independent  vari¬ 
ables.  The  single  intercorrelated  variable  chosen  for  regression  either  ex¬ 
plained  the  most  variation  in  dependent  (fish)  variables  or  was  the  most 
logical  relative  to  effects  documented  in  the  literature. 


Dependent  Variables 

State  conservation  agencies  provided  data  on  the  summer  catch  of  YOY 
fishes  in  a  variety  of  gears  including  seines  (in  Fort  Peck;  Oahe,  SD; 
Sharpe;  Frances  Case;  and  Lewis  and  Clark),  frame  nets  (Sakakawea  and 
Oahe,  ND),  and  gill  nets  (Sakakawea;  Oahe,  ND;  and  Francis  Case). 

Seines  were  100  by  9  ft1  with  0.25-in.  mesh,  and  frame  nets  were  3  by  4  ft 
with  0.25-in.  mesh  and  had  a  50-ft-long  lead.  Gill  nets  used  in  Sakakawea 
and  Oahe,  ND,  were  125  by  6  ft,  with  0.5-in.  monofilament  mesh.  Experi¬ 
mental  nets  used  in  Francis  Case  were  300  by  8  ft  with  six  50-ft  panels  of 
0.5-,  0.75-,  1.0-,  1.25-,  1.5-,  and  2.0-in.  mesh.  Most  sampling  was  in  Au¬ 
gust  or  September.  Data  were  transformed  by  taking  the  base  10  loga¬ 
rithm  of  one  plus  catch  and  averaged  to  obtain  one  value  per  reservoir, 
gear,  species,  and  year.  We  assigned  a  catch  of  zero  for  species  missing 
from  all  samples  in  a  year  if  it  was  captured  in  other  years.  The  YOY 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  given  on  page  vi. 
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LOG(ADULT  WALLEYE  CATCH) 


catch  of  walleye  in  Lake  Sakakawea  was  adjusted  so  that  only  data  on  nat¬ 
urally  produced  fish  were  used  for  regression  analysis.  The  catch  of  all 
YOY  walleye  was  multiplied  by  one  minus  the  fraction  that  was  stocked 
(number  marked/total  YOY  catch).  This  adjustment  was  possible  because 
stocked  walleye  were  intensively  marked  by  the  North  Dakota  Game  and 
Fish  Department  in  a  stocking-evaluation  project.  Catch  statistics  for 
YOY  fish  in  all  six  reservoirs  are  presented  in  Table  9. 

The  catch  of  YOY  fish  in  summer  was  used  as  a  dependent  variable  rep¬ 
resenting  reproductive  success  for  two  reasons.  First,  State  resource  agen¬ 
cies  for  Montana,  North  Dakota,  and  South  Dakota  indicated  that  the 
relative  abundance  of  YOY  fishes  is  a  fairly  reliable  indicator  of  future 
year-class  strength  in  these  reservoirs  (e.g..  Figure  1).  Although  an  abun¬ 
dant  cohort  of  YOY  fish  may  not  always  survive  to  create  a  strong  year 
class  (Fourt  1978),  the  probability  is  much  higher  than  when  few  VOY  are 
produced.  Second,  YOY  catch  and  water-level  changes  that  potentially  af¬ 
fect  catch  are  measured  in  the  same  year.  By  contrast,  the  catch  of  2-year- 
old  and  older  fish,  which  should  more  accurately  reflect  the  density  of 
harvestable  fish,  would  have  to  be  lagged  2  to  7  years  to  match  them  with 
hydrologic  conditions  that  may  have  produced  them.  Without  accurate 


Figure  1 .  Adult  walleye  catch  as  a  function  of  YOY  walleye  catch  in  earlier  years  in  Lake 
Sakakawea  (r2  =  0.71,  P  =>  0.0006) 
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age  and  growth  data,  this  lagging  process  could  have  an  error  as  large  or 
larger  than  the  one  potentially  resulting  from  the  assumption  that  densities 
of  harvestable  fish  are  proportional  to  YOY  densities  in  earlier  years.  The 
catch  of  older  fish  also  may  be  affected  by  factors  such  as  density-depen- 
dent  growth,  natural  mortality,  and  fishing  mortality  that  can  obscure  first- 
year  effects  of  water  levels. 


Correlation  and  Regression  Analyses 

Correlation  and  multiple-regression  analyses  were  used  to  find  the  best 
combinations  of  hydrologic  variables  for  predicting  YOY  catch  from  his¬ 
toric  data.  Dependent  catch  variables  were  matched  by  reservoir  and  year 
with  hydrologic,  weather,  and  covariate  stocking  variables.  We  generated 
Pearson  product-moment  correlation  matrices  with  the  CORR  Procedure 
and  regression  models  with  the  REG  Procedure  (MAXR  option)  of  the  Sta¬ 
tistical  Analysis  System  (SAS  Institute  Inc.  1988a  and  b).  Regression 
equations  were  evaluated  based  upon  statistics  such  as  the  significance 
level  of  the  model  and  parameter  estimates,  the  change  in  mean  square 
error  as  new  variables  were  added,  and  the  coefficient  of  determination  (r2 
or  R2).  We  did  not  accept  multiple-regression  models  if  the  relation  of 
YOY  catch  to  an  independent  variable  differed  (positive  or  negative)  from 
what  was  observed  in  correlation  analysis.  Most  importantly,  equations 
had  to  be  biologically  realistic  compared  with  known  ecological 
mechanisms. 

Results  of  two  rounds  of  correlation  and  regression  analyses  were  pre¬ 
sented  to  the  Reservoir  Fisheries  Task  Group.  After  the  first  round  of  anal¬ 
yses,  we  selected  regression  equations  meeting  the  criteria  described  in 
the  previous  paragraph.  Next,  we  picked  several  indicator  species  for 
each  reservoir  from  among  the  species  with  significant  relations  to  hydro- 
logic  variables.  Indicator  species  were  selected  because  they  were  import¬ 
ant  to  the  fishery,  were  highly  affected  by  operations,  or  represented  a 
distinct  spawning  strategy  (e.g.  riverine,  pelagic,  nesting,  or  broadcast  on 
vegetation).  Significant  relations  were  not  obtained  for  fish  representing 
every  spawning  strategy.  A  list  of  equations  was  presented  to  all  data  con¬ 
tributors  for  concurrence  on  the  best  equations  to  use  in  predictive 
software. 
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Postprocessing 


Background 

We  had  to  reduce  multi-reservoir,  -gear,  and  -species  predictions  into 
indicators  of  impacts  of  operational  alternatives  on  fish  reproduction  for 
the  entire  system  of  reservoirs  for  the  period  of  record.  The  MRD  needed 
this  data-reduction  process  to  develop  a  method  to  quickly  evaluate  hun¬ 
dreds  of  possible  operating  alternatives.  We  had  to  combine  predictions 
for  six  reservoirs,  three  sampling  gears,  and  from  two  to  five  species  per 
reservoir.  The  researchers  made  and  compared  predictions  (by  reservoir, 
gear,  and  species)  for  four  distinctly  different  operational  alternatives. 

The  integrated  model  uses  hydrologic  output  of  MRD’s  LRS  model  to 
calculate  annual  values  of  RI  (reproductive  index  from  1900  to  1990). 

The  RI  is  considered  the  best  available  index  for  impact  assessment  be¬ 
cause  models  were  based  upon  empirical  catch  data  for  several  species  of 
fish. 


Reproductive  index 

The  RI  was  calculated  in  a  five-step  process.  First,  we  made  predic¬ 
tions  of  YOY  catch  by  reservoir,  gear,  species,  and  year  (1900-1990).  Sec¬ 
ond,  we  standardized  predictions  by  dividing  predicted  catch  by  the 
maximum  observed  catch  for  the  same  species,  gear,  and  reservoir,  despite 
the  alternative,  so  each  species  was  weighted  equally.  Third,  standardized 
predictions  were  weighted  by  gear-specific  factors  and  area  of  habitat  to 
produce  an  index  to  the  total  number  of  YOY  fish  of  each  species  by  reser¬ 
voir,  gear,  and  year.  Standardized  seine  catches  were  divided  by  0.073, 
i.e.,  hectares  sampled  in  a  single  quarter  or  90-deg-arc  haul  and  multiplied 
by  the  mean  area  overlying  depths  of  0  to  30  ft  (assumed  depths  of  YOY 
habitat)  for  1  to  2  months  of  summer.  Months  included  June  and  July 
(Francis  Case  and  Lewis  and  Clark  Lakes),  July  (Lake  Sharpe),  July  and 
August  (Lake  Oahe),  and  August  and  September  (Sakakawea  and  Fort 
Peck  Lakes).  Standardized  gill-  and  frame-net  catches  were  divided  by 
0.1  (assuming  each  net  sampled  0.1  ha)  and  multiplied  by  area  over  0-  to 
30-ft  depths  in  the  months  listed  above.  Fourth,  we  summed  standardized 
indices  by  reservoir  and  year,  combining  different  gears  and  species. 

Fifth,  we  standardized  the  93  annual  RI  values  by  reservoir  by  dividing 
each  by  the  largest  RI  observed  for  that  reservoi”  "nder  any  alternative,  so 
each  reservoir  was  weighted  equally.  These  standardized  indices  were 
summed  by  year  to  index  fish  reproduction  for  the  entire  six-reservoir  sys¬ 
tem.  The  weighting  of  YOY  catch  by  habitat  area  in  each  reservoir  in  the 
third  step  and  differences  in  the  number  of  species  per  reservoir  forced  us 
to  standardize  by  reservoir  a  second  time  in  the  fifth  step. 
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We  weighted  predictions  by  reservoir  area  overlying  depths  of  0  to 
30  ft  to  account  for  differences  in  resource  size  among  alternatives  that  re¬ 
sulted  from  different  pool  levels.  Predicted  geometric  mean  catch  is  an  in¬ 
dicator  of  fish  density,  not  of  the  total  number  of  fish.  Total  number  is 
related  to  density  and  surface  area.  For  example,  a  reservoir  with  low- 
pool  elevations  during  drought  may  have  the  same  density  (number  per 
unit  area)  of  YOY  fish  as  it  does  when  the  basin  is  full.  However,  the 
total  number  of  YOY  present  would  be  higher  at  full  pool  because  there  is 
more  area  supporting  YOY  fish. 

Our  assumption  that  gill  and  frame  nets  sample  0. 1  ha  may  not  be  accu¬ 
rate,  but  0.1  is  a  constant  applied  to  all  gill-  and  frame-net  predictions,  de¬ 
spite  the  operational  alternative  evaluated.  Therefore,  the  procedure  is  no 
different  from  weighting  catch  by  surface  area  overlying  depths  of  0  to 
30  ft.  The  area  sampled  by  a  passive  gear  varies  greatly  because  of  fac¬ 
tors  affecting  fish  activity  and  movement.  Our  use  of  a  constant  sample 
area  was  more  to  show  that  the  quantity  was  unknown  than  to  assign  an  av¬ 
erage  value. 
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3  Results 


Area-  and  Volume-Elevation  Relations 

Quadratic  area-  and  volume-elevation  relations  (Table  1)  were  useful 
for  estimating  area  and  volume  from  elevation  data  recorded  from  1976  to 
1988.  Coefficients  of  determination  (R2)  for  all  quadratic  equations  ex¬ 
ceeded  0.99.  These  equations  would  be  less  accurate  for  data  collected  be¬ 
fore  1976  or  after  1988  because  of  sedimentation,  which  affects  area  and 
capacity  at  different  elevations.  Quadratic  area-elevation  equations  were 
better  predictors  of  surface  area  than  the  first  derivative  of  the  elevation- 
volume  equations,  which  are  sometimes  used  to  derive  area  estimates. 


Effects  of  Stocking 

Small  sample  sizes  and  poor  correlations  between  YOY  fish  catch  and 
most  stocking  variables  kept  us  from  using  stocking  variables  as  cov¬ 
ariates  in  regression  analyses.  In  this  study,  stocking  variables  limited  the 
number  of  years  that  could  be  included  in  a  model  and  reduced  most  sam¬ 
ple  sizes  to  five  or  less.  Stocking  records  seldom  provided  more  than 
5  years  of  data  for  any  species.  The  sample  size  of  a  multiple-regression 
model  is  determined  by  the  number  of  observations  in  which  every  inde¬ 
pendent  variable  has  a  nonmissing  value.  Observations  that  include  a 
missing  vnlue  for  any  independent  variable  in  a  model  are  dropped  from 
the  analysis. 

Correlations  shcw-d  that  fingerling  walleye  stocking  is  a  legitimate 
covariate,  although  sample  sizes  were  small.  We  found  positive  correla¬ 
tions  of  YOY  walleye  catch  with  stocking  variables  for  walleye  finger- 
lings  (but  not  fry)  in  Fort  Peck  Lake  (r  =  0.54,  P  =  0.1348,  N  =  8),  Lake 
Sakakawea  (r  =  0.76,  P  =  0.0304,  N  =  8),  and  Lake  Oahe,  SD  (r  =  0.62, 

P  =  0.0953,  N  =  8).  Insufficient  data  were  available  to  look  for  correla¬ 
tions  between  YOY  walleye  catch  and  stocking  in  Lake  Sharpe,  Lake  Fran¬ 
cis  Case,  or  Lewis  and  Clark  Lake. 
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Results  suggest  that  predictions  of  YOY  catch  could  be  significantly  im¬ 
proved  by  accounting  for  stocking  variation,  either  by  using  stocking  as  a 
covariate  (sample  sizes  permitting)  or  by  adjusting  catch  data  when 
stocked  fish  were  marked.  Our  adjustment  of  YOY  walleye  catch  in  Lake 
Sakakawea  to  include  only  nonstocked  YOY  resulted  in  a  stronger  relation 
between  YOY  catch  and  change  in  area  from  April  through  June  than 
when  catch  consisted  of  both  stocked  and  naturally  produced  walleye. 
Change  in  area  from  April  through  June  was  the  most  important  determi¬ 
nant  in  both  cases,  but  eliminating  stocking  effects  increased  the 
equation’s  r2  from  0.25  to  0.58  and  reduced  its  probability  from  0.0295  to 
0.0002.  Most  years  of  sport-fish  stocking  by  resource  agencies  occurred 
during  the  drought  of  the  1980’s,  which  probably  increased  apparent  repro¬ 
ductive  success,  as  indicated  by  catches  of  YOY  sport  fishes. 


Effects  of  Weather 

Correlation  of  YOY  catch  with  weather  variables  yielded  few  consis¬ 
tent  or  useful  results,  and  weather  variables  were  not  included  in  regres¬ 
sion  analyses.  We  found  positive  correlations  of  storm  hours  from  April 
through  June  with  the  catch  of  YOY  white  bass  in  Lake  Oahe,  SD 
(r  =  0.47,  P  =  0.0365,  N  =  20),  and  Lake  Sharpe  (r  =  0.7826,  P  =  0.0001, 

N  =  18)  and  with  the  catch  of  YOY  yellow  perch  in  Lake  Francis  Case 
(r  =  0.58,  P  =  0.0467,  N  =  12).  By  contrast,  storm  hours  from  April 
though  June  were  negatively  correlated  with  catches  of  YOY  sauger 
(r  =  -0.52,  P  =  0.0285,  N  =  18)  and  gizzard  shad  (r  =  -0.41,  P  =  0.0939, 

N  =  18)  in  Lewis  and  Clark  Lake.  The  75th-percentile  wind  speed  from 
April  through  June  was  positively  correlated  with  YOY  catches  of  white 
bass  (r  =  0.58,  P  =  0.0996,  N  =  9),  white  crappie  (r  =  0.61,  P  =  0.0844, 

N  =  9),  and  yellow  perch  (r  =  0.43,  P  =  0.0575,  N  =  20)  in  Lake  Oahe,  SD, 
and  walleye  (r  =  0.55,  P  =  0.0968,  N  =  10)  in  Lake  Francis  Case,  and  giz¬ 
zard  shad  (r  =  0.58,  P  =  0.01 14,  N  =  18)  and  sauger  (r  =  0.41,  P  =  0.0891, 
N  =  1 8)  in  Lewis  and  Clark  Lake.  White  crappie  in  Lake  Sharpe  were  in¬ 
versely  correlated  with  the  75th-percentile  wind  speed  (r  =  -0.67, 

P  =  0.0476,  N  =  9).  The  only  two  correlations  of  mean  equilibrium  tem¬ 
perature  from  April  through  June  with  YOY  catches  had  opposite  trends, 
one  positive  (Lake  Oahe  yellow  perch  (r  =  -0.48,  P  =  0.0298,  N  =  20))  and 
the  other  negative  (Lake  Sharpe  white  crappie  (r  =  0.73,  P  =  0.0266, 

N  =  9)). 

We  thought  that  storm  hours  during  spawning  would  be  inversely  re¬ 
lated  to  reproductive  success  of  many  species  because  wind-induced  turbu¬ 
lence  could  disrupt  spawning,  strand  eggs  and  larvae  along  shorelines,  or 
increase  silt  deposition  and  mortality  of  nonpelagic  eggs.  Surprisingly, 
three  of  the  five  correlations  we  found  were  positive.  Localized  storm 
events  could  increase  nutrient  loadings  from  the  immediate  watershed  and 
thereby  increase  primary  and  secondary  production  and  therefore  YOY 
survival.  Catches  of  YOY  white  bass  were  positively  correlated  with  sub¬ 
basin  inflow  in  Lake  Oahe.  Adult  white  bass  spawn  in  tributaries. 
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However,  we  believe  more  species  would  be  affected  if  productivity  were 
the  underlying  cause.  Also,  3  significant  correlations  out  of  26  possible 
correlations  of  YOY  catch  with  storm  hours  are  not  much  above  the  level 
of  chance  (0.05).  The  two  negative  correlations  with  storm  hours  ob¬ 
tained  for  two  species  in  Lewis  and  Clark  Lake  seem  to  support  our  origi¬ 
nal  hypothesis  about  negative  effects.  However,  storms  are  less  likely  to 
be  a  problem  for  fish  in  a  narrow  impoundment  like  Lewis  and  Clark  Lake 
than  they  would  be  in  a  reservoir  with  a  large  “fetch,”  i.e.  distance  over 
which  wind  can  blow  uninterrupted  by  land  (e.g.  Fort  Peck,  Sakakawea, 
Oahe,  and  Francis  Case).  Negative  correlations  of  storm  hours  with  YOY 
gizzard  shad  and  sauger  catches  could  be  related  to  high  turbidity  intro¬ 
duced  to  Lewis  and  Clark  Lake  by  the  Niobrara  River  during  the  spawning 
season.  Catches  of  YOY  sauger  and  gizzard  shad  also  were  inversely  re¬ 
lated  to  subbasin  inflow. 

The  importance  of  weather  at  one  or  two  sites  near  a  reservoir  to  fish 
reproduction  throughout  the  same  reservoir  may  be  questionable,  because 
weather  can  be  highly  localized.  At  best,  data  from  such  weather  stations 
might  realistically  portray  effects  of  widespread  fronts.  However,  they 
also  would  record  local  episodic  events  that  did  not  occur  30  to  100  miles 
away  and  would  miss  similar  events  on  other  areas  of  a  lake,  especially 
large  lakes  like  Fort  Peck  Lake,  Lake  Sakakawea,  and  Lake  Oahe.  This 
might  explain  the  six  positive  correlations  of  the  frequency  of  winds  ex¬ 
ceeding  the  75th-percentile  wind  speed  with  catches  of  YOY  fishes.  The 
only  negative  relation  was  obtained  for  white  crappie  in  Lake  Sharpe,  a 
mainstream  reservoir  with  little  fetch.  Documented  effects  of  wind  on 
fish  reproduction  have  been  exclusively  negative  (Clady  and  Hutchinson 
1975,  Summerfelt  1975,  Clady  1976,  Aggus  1979).  The  best  and  perhaps 
only  way  to  document  effects  of  weather  would  be  to  continuously  moni¬ 
tor  weather  at  multiple  fish-sampling  sites. 


Daily  Versus  Monthly  Hydrologic  Data 

Correlation  and  regression  analyses  using  hydrologic  variables  derived 
from  daily  data  provided  little  or  no  improvement  in  predictive  capability 
over  variables  derived  from  monthly  data.  The  same  hydrologic  variables 
usually  were  significant  or  nonsignificant  despite  the  time-step,  probably 
because  all  data  were  reduced  to  one  number  per  year  to  match  with  fish¬ 
ery  data.  Maximum  change  and  the  CV  in  area  from  15  April  through 
15  May  were  variables  that  could  be  calculated  only  from  daily  dab;.  Our 
hypothesis  was  that  these  variables  would  explain  variation  in  YOY  catch 
because  of  negative  impacts  of  drops  in  water  level  during  spawning. 
However,  both  variables  were  positively  correlated  with  spring  increases 
in  area  and  with  the  catch  of  several  YOY  fishes  the  next  August.  Either 
these  variables  do  not  capture  effects  of  brief  (1  to  2  day)  episodic  drops 
in  water  level  during  spawning,  or  such  sporadic  events  do  not  affect  YOY 
fish  production  as  much  as  other  factors  that  occur  after  spawning  (Gas- 
saway  1970). 
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Regression  Analyses 


We  found  many  highly  significant  relations  by  regressing  the  geometric 
mean  catch  of  YOY  fishes  on  hydrologic  variables  derived  from  monthly 
data.  Equations  retained  for  development  of  predictive  software 
(Table  10)  survived  careful  scrutiny  to  eliminate  relations  that  could  not 
be  explained  by  mechanisms  documented  in  the  literature.  For  example, 
years  of  very  high  inflow  are  associated  with  greater  surface  area  absorb¬ 
ing  solar  insolation,  increased  inundation  of  terrestrial  areas,  high  nutrient 
loading  (Perrier,  Westerdahl,  and  Nix  1977;  Westerdahl  et  al.  1981),  and 
increased  primary  and  secondary  production  (Benson  and  Cowell  1967, 
Dussart  et  al.  1972,  Vollenweider  1975,  Ostrofsky  and  Duthie  1978, 
McCammon  and  von  Geldem  1979,  Grimard  and  Jones  1982).  When  vege¬ 
tation  in  the  fluctuation  zone  is  flooded,  some  fishes  are  afforded  opti¬ 
mum  spawning  and  nursery  habitat,  e.g.  yellow  perch  (Beckman  and  Elrod 
1971),  northern  pike  (Benson  1968,  Hassler  1970),  buffaloes  (Moen 
1974),  and  common  carp  (Gabel  1974),  that  enhance  their  survival  (Mar¬ 
tin  et  al.  1981).  The  literature  on  effects  of  water  levels  and  inundation  of 
vegetation  is  replete  with  references  to  above-average  reproduction  and 
the  development  of  strong  year  classes  of  fish  under  such  conditions  (Ben¬ 
son  1968;  Beckman  and  Elrod  1971;  Nelson  and  Walburg  1977;  Nelson 
1978;  Ploskey,  Aggus,  and  Nestler  1985;  Ploskey  1986).  Regression  equa¬ 
tions  for  indicator  species  in  Fort  Peck  Lake,  Lake  Sakakawea,  Lake 
Oahe,  and  Lake  Francis  Case  (Table  9)  are  typical  of  positive  responses  to 
above-average  inflow  and  water  levels.  Densities  of  these  YOY  fish  usu¬ 
ally  were  highest  in  high-water  years,  in  spite  of  substantial  dilution  by  in¬ 
creased  water  volume. 

Multivariable  models  sometimes  contain  seasonal  inflow  or  change-in¬ 
area  variables  with  negative  coefficients,  but  usually  other  hydrologic  vari¬ 
ables  in  the  model  had  more  effect  on  the  cumulative  response.  Effects  of 
change-in-area  variables  cannot  be  interpreted  solely  by  noting  signs  of  co¬ 
efficients.  Equations  with  positive  coefficients  but  typically  negative  val¬ 
ues  indicate  that  small  decreases  in  area  are  more  beneficial  than  large 
ones.  Conversely,  equations  with  negative  coefficients  for  variables  with 
typically  negative  values  indicate  that  large  decreases  in  area  would  be 
better  for  fish  reproduction  than  small  decreases. 

We  know  that  YOY  fish  can  be  physically  concentrated  by  greatly  re¬ 
duced  water  levels  or  flushed  from  run-of-river  impoundments.  These 
mechanisms  can  obscure  or  override  our  ability  to  see  true  increases  in 
YOY  fish  densities  that  might  result  from  increased  system  productivity. 
Equations  for  fish  in  the  two  run-of-the-river  reservoirs  (Sharpe  and 
Lewis  and  Clark  Lakes)  probably  reflect  greater  physical  flushing  of  YOY 
fish  in  wet  years  (Walburg  1971),  as  relations  of  YOY  catch  to  inflow  and 
flushing  rate  variables  were  consistently  negative. 
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In  the  three  largest  reservoirs,  predicted  catches  were  positively  related 
to  flushing  rate  at  normal-pool  elevations  but  inversely  related  to  it  at  low- 
pool  elevations.  Flushing  rate  (discharge/volume)  increases  greatly  when 
reservoir  volume  becomes  very  low  because  it  is  a  ratio.  We  reran  regres¬ 
sions  after  excluding  seasonal  flushing  rate  variables  for  Fort  Peck  Lake, 
Lake  Sakakawea,  Lake  Oahe,  and  Lake  Francis  Case.  We  substituted  sea¬ 
sonal  inflow  for  flushing  rate.  We  observed  consistent  positive  relations 
between  predicted  catches  and  inflow  at  all  pool  elevations.  We  also 
dropped  equations  based  solely  on  inflow  variables  for  Fort  Peck  Lake  be¬ 
cause  inflow  to  Fork  Peck  does  not  vary  with  system  operations  nor 
among  alternatives. 


Integrated  Model  Application  to  Operating 
Alternatives 

Four  system-operating  alternatives  explored  in  this  study  differed 
mainly  in  system  storage  for  the  four  largest  reservoirs  and  inflows  to  the 
two  run-of-the-river  reservoirs  during  drought  (Figure  2).  Two  alterna¬ 
tives  differed  very  little  in  seasonal  water-level  or  hydrologic  patterns  in 
most  years,  so  impacts  to  fish  reproduction  were  most  obvious  in  drought 
years  and  predicted  indices  rarely  were  >3.0  (Figure  3).  A  second  pair  of 
operating  alternatives  allowed  for  significant  variation  in  seasonal  hydro- 
logic  patterns  in  many  years.  These  alternatives,  which  provided  a  year  of 
high  water  to  one  of  the  three  largest  reservoirs  on  a  rotating  basis,  pro¬ 
duced  similar  reproductive  indices  in  most  years.  However,  the  alterna¬ 
tive  allowing  the  greatest  summer  drawdown  produced  six  exceptionally 
high  RI  values  (>  3.0;  Figure  4).  It  also  yielded  more  years  with  above-av¬ 
erage  indices  (19  years  with  indices  >  2)  than  the  alternative  which  lim¬ 
ited  drawdown  (13  indices  >  2).  These  results  are  significant  because  a 
strong  year  class  of  fish  may  persist  for  about  5  to  8  years,  and  a  strong 
year  class  of  sport  fish  may  dominate  catches  of  anglers  for  3  to  3  years. 
However,  the  limited-drawdown  alternative  had  7  years  with  average  indi¬ 
ces  that  exceeded  indices  of  the  large-drawdown  alternative  from  1930  to 
1945,  a  period  of  drought.  Indices  for  both  alternatives  were  similar  in  6 
of  the  15  drought  years,  and  the  large-drawdown  alternative  produced 
higher  indices  than  the  limited-drawdown  alternative  in  two  of  these  years. 

The  exceptionally  high  RI  predicted  in  some  years  for  the  large- 
drawdown,  environmental  alternative  (Figure  4)  resulted  from  wet  years 
coinciding  with  refill  that  followed  a  drawdown  year  in  Lake  Oahe. 
Changes  in  pool  elevations  that  yielded  significant  differences  in  the  RI 
under  large-drawdown  and  limited-drawdown  alternatives  took  place  over 
2  years.  Exceptionally  high  indices  predicted  for  1914,  1929,  1971,  and 
1986  under  the  large-drawdown  alternative  were  not  predicted  for  the  lim¬ 
ited-drawdown  alternative  (Figure  4),  because  the  extent  of  drawdown  in 
the  previous  year  was  much  less  under  the  limited  drawdown  alternative 
(Figures  :>  and  6). 
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Figure  2.  System  storage  in  all  Missouri  River  reservoirs  predicted  by  MRD's  LRS  Model  under  two  operational 
alternatives 
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Figure  4.  Predicted  indices  of  fish  reproduction  for  six  Missouri  River  reservoirs  under  large-seasonal-drawdown  (triangles  and 
solid  line)  and  limited-seasonal-drawdown  (diamonds  and  dashes)  alternatives 
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Figure  5.  Lake  Oahe  water  elevations  under  large-  and  limited-drawdown  operational 

alternatives.  The  large-drawdown  alternative  produi  sd  exceptionally  high  fish 
reproductive  indices  in  1914  and  1929  in  contrast  to  the  average  indices  gener¬ 
ated  by  the  limited-drawdown  alternative 
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ELEVATION  (mean  sea  level) 


Figure  6.  Lake  Oahe  water  elevations  under  large-  and  limited-drawdown  operational 

alternatives.  The  large-drawdown  alternative  produced  exceptionally  high  fish 
reproductive  indices  in  1971  and  1986  in  contrast  to  the  average  indices  gener¬ 
ated  by  the  limited-drawdown  alternative 
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Using  Individual  Regression  Models 


Single  equations  (Table  10)  can  be  used  to  forecast  the  YOY  catch  of 
select  species  of  fish  from  hydrologic  variables  (Table  2)  derived  from 
monthly  elevation,  inflow,  and  release  data.  End-of-month  data  can  be  ob¬ 
tained  from  MRD  or  U.S.  Geological  Survey  Surface  Water  Records. 

Three  steps  are  required  to  make  predictions.  First,  obtain  or  create  end- 
of-month  data  on  elevation,  inflow,  and  release  for  the  year  to  be  pre¬ 
dicted  and  the  previous  year  for  a  specific  reservoir.  Second,  select 
equations  for  the  reservoir  and  species  of  interest  (Table  10);  list  indepen¬ 
dent  variables  used  in  the  equations,  and  obtain  definitions  of  independent 
variables  (Table  2).  Third,  calculate  values  of  independent  variables  ac¬ 
cording  to  the  definitions;  substitute  calculated  values  in  the  appropriate 
regression  equation,  and  solve  the  equation  for  YOY  catch.  See  Chapter 
4,  Conclusions  and  Limitations,  for  information  on  predictions  using  val¬ 
ues  of  independent  variables  that  are  outside  the  range  of  empirical  values 
used  to  derive  the  equations. 

Regression  equations  in  Table  10  predict  the  base  10  logarithm  of  YOY 
catch  +  1,  which  can  be  converted  to  the  geometric  or  arithmetic  mean 
(Ricker  1975).  The  geometric  mean  ( 1  o,0*<c*teh)  -  1)  can  be  converted  to 
an  approximate  arithmetic  mean  using  the  following  formula: 

log,0(AM)  =  1.1518s2  (N  -  1)/N  +  log,0(GM) 

where  AM  is  the  arithmetic  mean,  GM  is  the  geometric  mean,  H  is  sample 
size,  and  £  is  the  standard  deviation  of  the  normally  distributed  base  10 
logarithms  of  catch  (Table  9). 
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4  Conclusions  and 
Limitations 


Models  can  be  used  to  evaluate  seasonal  operating  alternatives  or  differ¬ 
ences  in  system-operating  alternatives  (after  postprocessing  of  predic¬ 
tions),  as  described  in  previous  sections. 

Our  evaluation  of  four  system-operating  alternatives  suggests  that  pro¬ 
vision  of  a  year  of  high  water  to  one  of  the  three  largest  reservoirs  on  a  ro¬ 
tating  basis  yields  the  greatest  benefit  to  natural  fish  reproduction  for  the 
system.  Alternatives  that  limit  annual  drawdown  are  desirable  only  for  se¬ 
vere  drought  periods  when  the  fish  reproduction  and  reservoir  fisheries 
are  both  adversely  affected  by  low  water.  In  normal  water  years,  large 
summer  drawdown  followed  by  a  year  of  above-average  water  levels  can 
greatly  increase  fish  reproductive  success. 

Software  developed  in  this  study  allows  users  to  make  annual  or  multi¬ 
year  predictions  quickly,  but  the  present  version  does  not  screen  values  of 
independent  variables  to  make  certain  they  are  within  the  range  of  empiri¬ 
cal  data  used  to  derive  the  regression  equations.  Extrapolation  is  a  con¬ 
cern  primarily  for  users  making  predictions  of  YOY  catch  with  individual 
equations.  In  these  cases,  input  data  should  be  screened,  or  users  must  as¬ 
sume  that  relations  are  consistent  over  a  wider  range  of  values  of  indepen¬ 
dent  variables  than  ever  observed.  Concern  over  extrapolation  should 
partly  depend  on  how  far  a  value  is  out  of  range.  For  example,  values 
over  100  percent  out  of  range  might  be  considered  risky,  whereas  those  1, 
10,  20,  or  even  30  percent  out  of  range  may  be  believable.  Users  making 
predictions  with  equations  in  Table  10  should  compare  calculated  values 
of  independent  variables  to  the  maxima  and  minima  listed  in  Tables  3  to  8 
to  identify  years  when  predictions  may  be  suspect.  Users  could  assign  the 
minimum  or  maximum  (Tables  3  to  8)  to  an  outlying  value  to  assure  that 
predictions  are  within  the  original  range  of  values. 

Extrapolation  beyond  the  original  data  is  not  a  serious  problem  for  the 
integrated  model  that  uses  93  years  of  predicted  hydrology,  because  YOY 
predictions  are  standardized  by  reservoir  and  species  to  values  between 
zero  and  one.  Consequently,  a  prediction  from  a  single  equation  cannot 
overly  bias  the  composite,  annual  estimate  of  the  RI.  In  addition,  the 
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model  is  used  solely  to  compare  alternatives,  not  to  make  quantitative  pre¬ 
dictions.  Users  evaluating  system  alternatives  with  93  years  of  predicted 
hydrology,  including  the  extreme  drought  of  the  1930’s  and  the  wettest 
years  recorded,  likely  will  be  beyond  the  range  of  the  original  data  in 
some  years.  The  environmental  alternative  that  sought  to  provide  high 
pools  in  one  of  the  three  largest  reservoirs  every  third  year  and  allowed  a 
large  seasonal  drawdown  had  the  most  outliers  (Table  11).  Variables  in¬ 
dexing  inflow,  flushing  rate,  or  change  in  area  in  spring  were  the  most 
common  offenders. 
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Table  1 

Coefficients  in  quadratic  relations  between  surface  area  (acres)  or  volume 
(acre-ft)  and  elevation  (ft,  mean  sea  level)  in  six  mainstream  Missouri  River 
Reservoirs1 


Lake 

AO 

A1 

A2 

VO 

VI 

V2 

FP 

17443774.5664 

-17355.6704 

4.3156 

2522694031.9262 

•2433787.8962 

586.9796 

SA 

27276713.1370 

-32890.2021 

9.9136 

2967074723.5567 

•3507845.2696 

1029.7744 

OA 

23067376.9553 

-31945.4196 

11.0682 

-2710709.6077 

9252710 

SH 

10478680.7020 

-15883.9167 

6.0172 

725648789.2619 

-1071504.7706 

3952526 

FC 

4571151.6347 

-7984.9174 

3.4485 

653173916.9068 

-10363422346 

411.1322 

LC 

17935605.2458 

-30653.2061 

13.2697 

637158493.1287 

-1083328.8038 

460.4746 

Equations  have  the  farm:  AREA  -A0  +  A1  (ELEV)  +  A2(ELEV)2  Of  VOL  -  VO  +  VI  (ELEV)  +  V2(ELEV)2,  when*  AO, 
A1,  A2,  VO,  VI,  and  V2  are  coefficients  tabled  above,  ELEV  is  aiavation,  and  VOL  is  volume.  Lake  abbreviations 
are  as  Mows:  FP  -  Fort  Peck;  SA  -  Sakata wea;  OA  -  Oahe;  SH  -  Sharpe;  FC  -  Francis  Case;  and  LC  -  Lewis 
and  Clark. 
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Table  2 

Independent  hydrologic  variables,  definitions,  and  transformations 


Variable 

Definition 

Transformation 

Mean  discharge  at  an  upstream  gage  (inflow)  per  year  in  cms 

-iogto 

CASUSP 

Change  in  mean  area  (±  ha)  summer  to  spring,  i.e.,  mean  of  areas  at  the 
end  of  Apr,  May,  and  Jun  minus  the  mean  of  areas  at  the  end  of  Jut,  Jul, 
Aug,  and  Sep  in  the  previous  year 

-none 

CASUSP2 

Sum  of  changes  in  mean  area  (±  ha)  from  summer  to  spring  tor  two 
consecutive  years  (see  CASUSP  above) 

-none 

CASUSU 

Change  in  mean  area  (±  ha)  from  summer  to  summer,  i.e.,  mean  of  areas 
at  the  end  of  Jut,  Jul,  Aug,  and  Sep  minus  the  mean  of  areas  at  the  end 
of  Jun,  Jul,  Aug,  &  Sep  in  the  previous  year 

•none 

CASUSU2 

Sum  of  changes  in  mean  area  (±  ha)  from  summer  to  summer  for  two 
consecutive  years 

•none 

INFV4_5 

Total  inflow  volume  (mMons  of  cubic  meters)  from  31  Mar  through  31  May 

-toflio 

INFV46 

Total  inflow  volume  (mittons  of  cubic  meters)  from  31  Mar  through  30  Jun 

-logio 

XSBINF4_6 

Mean  of  subbasin  inflows  on  31  Mar,  30  Apr,  31  May,  and  30  Jun  (cms) 

-  logic 

XPA4  5 

Area  (ha)  associated  with  the  mean  of  elevations  on  31  Mar,  30  Apr,  and 

31  May  minus  area  at  an  elevation  30  ft  beiow  the  mean 

-logio 

XPA46 

Area  (ha)  associated  with  the  mean  of  aievalions  on  31  Mar,  30  Apr,  31 
May.  and  30  Jun  minus  area  at  an  elevation  30  ft  below  the  mean 

-log,o 

X20V4_6 

Volume  (milions  of  cubic  meters)  associated  with  the  mean  of  elevations 
on  31  Mar,  30  Apr,  31  May,  and  30  Jun  minus  volume  at  an  elevation  30  ft 
below  the  mean 

-•ogio 

CA4_S 

Change  in  area  (±  ha):  31  Mark)  31  May 

-none 

CA46 

Change  in  area  (±  ha):  31  Mar  to  30  Jun 

-none 

FR4_5 

Flushing  rate  from  Apr  through  May,  where  flushing  rale  is  the  total 
release  divided  by  mean  volume 

-logio 

FR4_6 

Flushing  rate  from  Apr  through  Jun,  where  flushing  rate  is  the  total 
released  divided  by  mean  volume 

-logio 

INFV6_9 

Total  inflow  volume  (millions  of  cubic  meters)  from  31  May  through  30  Sep 

-logio 

XPA7_11 

Area  (ha)  associated  with  the  mean  of  elevations  on  31  May,  30  Jun,  31 

Jul,  31  Aug,  30  Sep,  31  Oct,  and  30  Nov  minus  area  at  an  elevation  30  ft 
below  the  mean 

-logio 

XPA50_6_9 

*  -*w  (ha)  associated  with  the  mean  of  elevations  on  31  May,  30  Jun,  31 

Jul,  31  Aug,  and  30  Sep  minus  area  at  an  elevation  50  ft  below  the  mean 

-logio 

CA6_9 

Change  in  area  (±  ha):  31  May  through  30  Sep 

-none 

CA7_11 

Change  in  area  (±  ha):  30  Jun  through  30  Nov 

-none 

FR6_9 

Rushing  rate  from  Jun  through  Sep,  where  flushing  rate  is  the  total 
release  divided  by  mean  volume 

-logio 

Table  4 

Sample  sizes  (N),  minima,  lower  qualities,  medians,  means,  upper  quartiles,  and 
maxima  of  hydrologic  variables  derived  for  Lake  Sakakawea,  North  Dakota 


MG_WY  24 


CASUSP 


CASUSP2 


CASUSU  23 


CASUSU2  22 


1NFV4_6  24 


XSBINF4_6  24 


XPA4_6  24 


XPV4_6  24 


CA4_6  24 


INFV6_9  24 


XPA6_9  24 


XPA50J8J9  24 


CA6  9  24 


Minimum 


387.00 


-15103.00 


-24738.00 


-17808.00 


-24648.00 


5581.00 


324.70 


24234.00 


55099.94 


•1254.00 


5529.00 


24446.98 


55099.94 


-5852.00 


Quart!  I* 


560.85 


-10474.00 


-15151.00 


-6840.00 


-10224.00 


7309.40 


462.35 


26668.94 


62345.20 


3092.00 


8442.70 


27345.20 


6234520 


785.00 


Median 

ai - 

wan 

Ud Mr 

Quart!  la 

725.93 

674.54 

820.02 

-6107.00 

•638826 

-1941.00 

-13920.00 

-12692.50 

•6292.00 

-1090.00 

-1347.30 

3934.00 

-59.50 

-2593.05 

4200.00 

9673.96 

9553.90 

12800.03 

635.00 

631.48 

870.09 

27352.95 

27079.71 

27801.00 

64279.46 

63139.35 

65005.98 

6327.00 

7474.88 

11263.50 

1112128 

10726.49 

13542.19 

2811920 

27665.49 

28409.51 

64279.46 

63139.35 

65005.98 

583320 

5805.08 

11001.50 

Maximum 


101820 


189.00 


-585.00 


13318. 


14395.00 


14778.99 


1175.50 


28457.98 


.94 


26694.00 


20395.99 


29049.99 


66806.94 


15664.00 


Note:  Variable  abbreviations  are  defined  in  Table  2. 


Table  5 

Sample  sizes  (N),  minima,  lower  quartiies,  medians,  means,  upper  quartlles,  and 
maxima  of  hydrologic  variables  derived  for  Lake  Oahe,  North  and  South  Dakota 


pwei 

L ower 
Quartiia 

Median 

»i- - 

Mean 

E 

Maximum 

27 

512.50 

620.50 

713.50 

738.19 

881.10 

1075.30 

CASUSP 

26 

-12545.00 

-7087.00 

ICS3 

1790.73 

11930.00 

21845.00 

CASUSP2 

25 

-22136.00 

-10876.00 

2176.00 

2784.16 

10582.00 

34994.00 

CASUSU 

26 

-17350.00 

-9778.00 

26150 

1139.08 

12153.00 

26204.00 

CASUSU2 

25 

-33614.00 

-12760.00 

1344.88 

12621.00 

31305.00 

INFV4_6 

27 

5950.00 

6768.01 

80C5.00 

8663.69 

10872.99 

15041.01 

XSBINF4_6 

27 

0.00 

42-50 

143.00 

7952 

313.70 

505.80 

XPA4_6 

27 

19957.00 

25895.03 

27812.97 

2702159 

28552.02 

28897.01 

CA4_6 

27 

-10412.00 

-2362.00 

4543.82 

10172.00 

18515.00 

INFV6.9 

27 

6357.00 

7355.00 

8953.00 

9480.00 

12065.99 

18373.99 

XPA6J9 

27 

19598.00 

26037.98 

27489.00 

26909.51 

28796.98 

29047.99 

CA6_9 

27 

-14326.00 

•10805.00 

-6721.00 

-6184.74 

-3810.00 

9906.00 

XPA7_1 1 

27 

19327.00 

25702.01 

27210.03 

26596.35 

28208.02 

28687.97 

CA7_1 1 

27 

-20873.00 

•13246.00 

-10436.00 

-9767.41 

-7004.00 

14626.00 

1 


Variable  abbreviations  are  defined  in  Table  2. 


Table  6 

Sample  sizes  (N),  minima,  I  *  quartiles,  medians,  means,  upper  quartiles,  and 
maxima  of  hydrologic  variable*  derived  for  Lake  Sharpe,  South  Dakota 


Variable1  N  Minimum  Quartito  Mad  ten  linn  Quart!  le  Maximum 


MG_WY 

28 

456.90 

670.82 

691.80 

786.84 

1044.40 

INFV4_6 

27 

mm 

5920.99 

7048.00 

6862.63 

7866.99 

10725.01 

XSBINF4_6 

28 

0.00 

2.17 

14.89 

8.44 

21.05 

74.20 

LFFI4_6 

27 

1.91 

2.74 

3.30 

3.22 

3.73 

4.94 

INFV6_9 

28 

3439.00 

9411.41 

10241.08 

10388.84 

10388.84 

16265.01 

FR6_9 

28 

3.84 

4.41 

•*.d3 

— 

6.00 

6.00 

294.56 

1  Variable  abbreviations  are  defined  in  Table  2. 


Table  7 

Sample  sizes  (N),  minima,  lower  quartiles,  medians,  means,  upper  quartiies,  and 
maxima  of  hydrologic  variables  derived  for  Lake  Francis  Case,  South  Dakota 


VartaMa1 


MG_WY 

28 

CASUSP 

27 

CASUSP2 

26 

CASUSU 

27 

CASUSU2 

26 

INFV4  6 

| 

28 

XS8INF4_6 

28 

XPA4_5 

28 

CA4_5 

28 

INFV6_9 

28 

XPA6_9 

28 

CA6J9 

28 

.50 


-2485.00 


-4327. 


4921.00 


0.00 


Uppw 

Mtdian 

Maan 

Quartit* 

Maximum 

656.07 


-290.00 


691.55  722.13 


820.01  1082.90 


525.00  667.96  1657.00  5126.00 


-2568.00  -403.00  1080.50  1363.04  2483. 


-3689.00  -982.00 


-1005.00 


82.93  1483.00  3067.00 


238.65  1413.00  5538.00 


7579.49  7463.40  8578.48  11590.99 


10134.01  10441.49  10511.50  10518.01  10614.49  10791.00 


•2104.00  543.50 


867.64  216450  6554.00 


8374.00  9906.00  1048750  11424.71  12620.50  17104.00 


10005.00  1041050  1047250  10468.86  10575.00  10810.00 


•6192.00  -3880.00  -2529.00  -2671.14  -1483.50 


1 


Variable  abbreviations  are  defined  in  Table  2. 


Table  8 

Sample  sizes  (N),  minima,  lower  quart! les,  medians,  means,  upper  quart! les,  and 
maxima  of  hydrologic  variables  derived  for  Lewis  and  Clark  Lake,  South  Dakota 


Variable1  N  Mini  mom  Quart!  la  Mad  Ian  Maan  Quartlla  Maximum 


MG_WY 

28 

580.30 

711.68 

770.01 

779.24 

870.72 

1121.20 

INFV4_6 

28 

5586.99 

7418.51 

8018.53 

7933.82 

8800.38 

11218.00 

XSBINF4_6 

28 

38.60 

67.10 

81.40 

85.78 

113.89 

195.70 

FR4_6 

28 

11.54 

15.28 

17.00 

17.05 

18.73 

23.83 

INFV6_9 

28 

10070.01 

11280.95 

12192.49 

12503.92 

13703.94 

18325.98 

FR6_9 

28 

18.82 

21.54 

23.27 

24.46 

26.74 

35.98 

1  Variable  abbreviations  are  defined  in  Table  2. 


Table  9 


Distribution  statistics  for  log-transformed  catches  of  young-of-year  fishes  by 
reservoir,  species,  and  gear 


Fort  Pack  northern 


pike 

in  seines;  Variable  * 

tofli 

number 

N 

Mean 

Std  Dev 

Skewness 

CV 

17 

0.178985 

0247881 

2.508406 

138.4925 

Sum 

Variance 

Kurtosis 

Std  Mean 

3.042749 

0.061445 

5205228 

0.06012 

Quantiles 

100%  Max 

0239519 

95% 

0.939519 

75%  03 

0.173186 

90% 

0.672098 

50%  Med 

0.079181 

10% 

0.041393 

25%  Q1 

0.041393 

5% 

0.041393 

0%  Min 

0.041393 

r II  n — 

Mn  KfCa 


number 


StdOev 

Skewness 

CV 


Quantiles 
100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


Sid  Dev 

Skewness 

CV 


Quantiles 
100%  Max 
75%  Q3 
50%  Med 
25%  Q1 
0%Min 


N 

Mean 
Std  Dev 
Skewness 
CV 


Quantiies 
100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


14 

0.141113 

0.062443 

-0.20841 

4425036 


0255273 

0.173186 

0.159657 

0.079181 

0.041393 


Fort  Peck  yellow 


17 

1.658385 

0.598484 

-0.74726 


2255643 

2.144574 

1.79588 

1.089905 

0.462398 


Sakakaweawal 


19 

0.040194 

0.031005 

1.148238 

77.13984 


0.112605 

0.66699 

0.026125 

0.017033 

0.008174 


Sum 
Variance 
Kurtosis 
Std  Mean 


In  seines;  Variables 

>Pfl< 

number 

Sum 
Variance 
Kurtosis 
Std  Mean 


1275576 


-0.46756 

0.016689 


0255273 

020412 

0.041393 

0.041393 


28.19255 

0258183 

-021334 

0.145154 


In  gill  nets;  Variable  s 

log( 

number 

Sum 
Variance 
Kurtosis 
Std  Mean 


0.763683 

0.000961 


0.007113 


0.112605 

0.101059 

0.009451 

0.008174 


(Shmt1ot7) 


1  Definitions  of  abbreviated  variables  are  as  follows:  N  -  sample  size;  Std  Dev  -  standard  deviation;  CV  -  coefficient  of 
variation;  Std  Mean  «  standard  error  of  the  mean;  Max  »  maximum;  03  »  75th  percentile;  Med  -  median;  Q1 «  25th 
percentile;  Min  -  minimum. 


Table  9  (Continued)  _ _ 


Sefcafcawee  crapple  In  frame  nett;  Variable  » log  (number  per  hour*1) 


Sid  Dm 

Skewness 

CV 


Quantiles 
100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


19 

0.095631 

0.069531 

1.266363 

93.62149 


0.31597 

0.149219 


0.021189 

0.004321 


Sum 
Variance 
Kurtosis 
Std  Mean 


1.816995 

0.008016 

0.96743 

0.02054 


0.31597 

0274158 

0.017033 

0.004321 


Oahe,  NO,  walleye  tn  gill  nets;  Variable  a  log  (number  per  hour  ♦  1) 


N 

Mean 

Std  Dev 

Skewness 

CV 

17 

0.026226 

0.041458 

2.888644 

158.0769 

Sum 

Variance 

Kurlosis 

Std  Mean 

0.445847 

0.001719 

9.056593 

0.010055 

Quantiles 

100%  Max 

0.168055 

95% 

0.168055 

75%  Q3 

0.02214 

90% 

0.073755 

50%  Med 

0.010342 

10% 

0.002986 

25%  Q1 

0.006124 

5% 

0.000738 

0%  Min 

0.000738 

N 

Mean 
Sid  Dev 
Skewness 
CV 


Quantiles 
100%  Max 
75%  Q3 
50%  Med 
25%  Q1 
0%  Min 


N 

Mean 
Std  Dev 
Skewness 
CV 


100%  Max 
75%  Q3 
50%  Med 
25%  Q1 
0%  Min 


Oahe,  ND,  white 


17 

0.127136 

0.135155 

1.178775 

106.3074 


0.466853 

0.225361 

0.072654 


I  nets;  Variable  « log  (number  per  hour  ♦  1) 


1 

Sum 

2.161315 

Variance 

0.018267 

Kurtosis 

0.866936 

Std  Mean 

0.03278 

95% 

0.466853 

90% 

0.313572 

10% 

0 

5% 

0 

Oahe,  ND,  crappie  In  small  frame  nett;  Variable  a  log  (number  per  hour  + 1) 


4.775437 

0.084429 

0.09985 

0.070473 


17 

0280908 

0290567 

1.148339 

103.4384 

Sum 

Variance 

Kurtosis 

Std  Mean 

0.871666 

95% 

0.312389 

90% 

0.183298 

10% 

0.03153 

5% 

0.003719 

0.871666 

0.841422 

0.01695 

0.003719 


(Sh*M2of7)\ 


Table  9  (Continued) 


Oaha,  NO,  yellow  parch  in  small  trams  nets;  Variable  =  log  (number  par  hour  ♦  1) 


N 

Mean 

StdDev 

Skewness 

CV 

17 

0.182641 

0261334 

1.903691 

143.0864 

Sum 

Variance 

Kurtoeis 

Sid  Mean 

3.104894 

0.068296 

2.766081 

0.063383 

Quantiles 

100%  Max 

0.815036 

95% 

0.815936 

75%  03 

0.234796 

90% 

0.815452 

50%  Mad 

0.042733 

10% 

0.004536 

25%  Q1 

0.032498 

5% 

0.000651 

0%  Min 

0.000651 

Oaha,  SO,  northern  pika  In  aeinea;  Variable  »  log  (number  par  haul  ♦  1) 


N 

Mew 

StdDev 

Skewness 

CV 

20 

0.027704 

0.054545 

2.541252 

196.8832 

Sum 

Variance 

Kurtosis 

Std  Mean 

0.554082 

0.002975 

6.08148 

0.012197 

Quantiles 

100%  Max 

020412 

95% 

0.175124 

75%  Q3 

0.020775 

90% 

0.112655 

50%  Med 

0 

10% 

0 

25%  Q1 

0 

5% 

0 

0%  Min 

0 

Oahe,  SO,  walleye  in  aaktaa;  Variable  *  log  (number  par  haul  + 1) 


N 

Mean 

StdDev 

Skewness 

CV 

20 

0.153098 

0.19212 

1.614917 

125.4878 

Sum 

Variance 

Kurtosis 

Std  Mean 

3.061968 

0.03691 

2.330558 

0.042959 

Quantiles 

100%  Max 

0.69897 

95% 

0.588046 

75%  Q3 

0.190106 

90% 

0.454243 

50%  Med 

0.096562 

10% 

0 

25%  Q1 

0 

5% 

0 

0%  Min 

0 

Oaha,  SO,  white  base  in  aeinea;  Variable  =  log  (number  par  haul  ♦  1) 


N 

Mean 

StdDev 

Skewness 

CV 

20 

1.344825 

0.422464 

0.337966 

31.41404 

Sum 

Variance 

Kurtosis 

Std  Mean 

26.89649 

0.178476 

-03869 

0.094466 

Quantiles 

100%  Max 

2206826 

95% 

2.099247 

75%  Q3 

1.607157 

90% 

1.98922 

50%  Med 

1.30533 

10% 

0.757672 

25%  Q1 

1.011832 

5% 

0.707487 

0%Min 

0.69897 

Table  9  (Continued) 


Oahe,  SO,  white  crappie  in 


;  Variable  *  log  (number  per  haul  ♦  1) 


StdDev 

Sinmms 

CM 


Quantiles 
100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


9 

0.262832 

0.472856 

2.48612 

179.9079 


1.462398 

0.30103 

0 

0 

0 


Sum 

Variance 

Kurtosis 

StdMean 


2.365488 

0223592 

6.650386 

0.157619 


Oahe,  SO,  yellow  perch  In  seines;  Variable  ■  log  (number  per  haul  + 1) 


N 

20 

Mean 

1.437807 

Sum 

28.75615 

StdDev 

0.716007 

Variance 

0212666 

Skewness 

0.099228 

Kurtosis 

CV 

49.79855 

StdMean 

0.160104 

Quantiles 

100%  Max 

2.76809 

95% 

2.630704 

75%  Q3 

1.925718 

90% 

226455 

50%  Med 

1.380091 

10% 

0.477121 

25%  Q1 

0.80103 

5% 

0.477121 

0%  Min 

0.477121 

Sharpe  gizzard  shad  In  seines;  Variable  =  log  (number  per  haul  +  1) 

N 

18 

Mean 

2.321466 

Sum 

41.78639 

StdDev 

0.459592 

Variance 

0211224 

Skewness 

0.008068 

Kurtosis 

>1.09756 

CV 

19.79747 

StdMean 

0.108327 

Quantiles 

100%  Max 

3.131939 

95% 

3.131939 

75%  Q3 

2.659916 

90% 

2.904174 

50%  Med 

2217543 

10% 

1.675778 

25%  Q1 

1.965202 

5% 

1.61595 

0%  Min 

1.61595 

Sharpe  freshwater  drum  In  seines;  Variable  =  log  (number  per  haul  ♦  1) 


N 

Mean 

StdDev 

Skewness 

CM 


100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


18 

0.420005 

0236 

0.351766 

79.99886 

Sum 

Variance 

Kurtosis 

StdMean 

7.560098 

0.112896 

-1.13338 

0.079196 

1.079181 

95% 

1.079181 

0.69897 

90% 

0.845098 

0.406457 

10% 

0.041393 

0.079181 

5% 

0 

0 

(Sh—t4ot7)\ 


Table  9  (Continued) 


Sharp*  walleye  In  —In—;  Variable  =  log  (number  per  haul  ♦  1) 


N 

Mean 

StdOev 

Skewness 

C  V 

18 

0.66360 

0.30619 

156178 

46.14085 

Sum 

Variance 

Kurtosis 

Std  Mean 

11.94478 

0.093752 

2.540681 

0.07217 

Quanta— 

100%  Max 

1531479 

95% 

1.531479 

75%  03 

0.845096 

90% 

0.954243 

50%  Med 

0.579633 

10% 

0.30103 

25%  01 

0.477121 

5% 

0578754 

0%  Min 

0.278754 

Sharpe  white  ba—  in  — Inaa;  Variable  *  log  (number  par  haul  + 1) 


N 

18 

Mean 

0.633824 

Sum 

11.40884 

Std  Dev 

0.502183 

Variance 

0552187 

Skewness 

0.579233 

Kurtosts 

-0.99596 

CV 

79.23059 

Std  Mean 

0.118366 

Quanta— 

100%  Max 

1.544068 

95% 

1.544068 

75%  03 

1.041393 

90% 

1.447158 

50%  Med 

0.560287 

10% 

0.079181 

25%  Q1 

0.146128 

5% 

0.079181 

0%  Min 

0.079181 

Francis  Ca—  gizzard  altad  in  aeln— ;  Variabla  *  log  (number  per  haul  ♦  1) 

N 

22 

Mean 

1.429021 

Sum 

31.43846 

Std  Dev 

0.70904 

Variance 

0502738 

Skewness 

0.007643 

Kurtosis 

-058758 

CV 

49.6172 

Std  Mean 

0.151168 

Quanta— 

100%  Max 

2531479 

95% 

2.454845 

75%  Q3 

1.907551 

90% 

2.424882 

50%  Med 

1531719 

10% 

0.795045 

25%  Q1 

0.981229 

5% 

057277 

0%  Min 

0.01536 

Francis  Case  walleye  In  gill  nets;  Variable  =  log  (number  per  net  night  *  1) 


N 

Mean 

Std  Dev 

Skewness 

CV 

10 

0.539595 

0.3296 

•0.11496 

61.08276 

Sum 

Variance 

Kurtosis 

Std  Mean 

5.395955 

0.108636 

-0.73228 

0.104229 

Quanta— 

100%  Max 

0.995635 

95% 

256*216 

75%  Q3 

0.748188 

90% 

2567172 

50%  Med 

0.504938 

10% 

0.931458 

25%  Q1 

0%  Min 

0.342423 

0 

5% 

0.863561 

(Sh—t5of7)\ 


Table  9  (Continued) 


Francis  Cm*  whit#  bass  In  seines;  Variable  *  k>fl  (number  per  haul  ♦  1) 


StdDsv 

Skewness 

CV 


Quantiles 
100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


21 

1.491915 

0.551972 

0.150358 

36.99754 


2.517196 

1.83089 

1.407136 

1.056486 

0.390228 


Sum 
Variance 
Kurtoeis 
Sid  Mean 


31.33021 

0.304673 

-0.50528 

0.12045 


2.369216 

2.267172 

0.931458 

0.863561 


Francis  Case  white  crappie  In  seines;  Variable  ■  log  (number  per  haul  +  1) 


StdOev 

Skewness 

CV 


Quantiles 
100%  Max 
75%  03 
50%  Med 
25%  Q1 
0%  Min 


0.491308 

1.755407 

152.1898 


1.63731 

0.477121 


Sum 
Variance 
Kurtosis 
Sid  Mean 


6.779338 

0.241383 

2.185588 

0.107212 


1.380211 

1.189041 

0 

0 


Francis  Case  yellow  parch  In 


;  Variable  *  log  (number  per  haul  +1) 


N 

22 

Mean 

1.08346 

Sum 

23.83612 

StdOev 

0.52817 

Variance 

0.27896 

Skewness 

0.12191 

Kurtosis 

-0.23661 

CV 

48.74848 

Std  Mean 

0.11260 

Quantiles 

100%  Max 

2.05490 

95% 

2.032728 

75%  03 

1.38806 

90% 

1.748188 

50%  Med 

1.03984 

10% 

0.519566 

25%  Q1 

0.69897 

5% 

0.374235 

0%  Min 

0 

Lewis  and  Clark  gizzard  shad  in  seines;  Variable  =  log  (number  par  haul  ♦  1) 

N 

18 

Mean 

1.968775 

Sum 

35.43794 

Std  Dev 

0.613034 

Variance 

0.375811 

Skewness 

0.025618 

Kurtosis 

-0.89553 

CV 

31.13786 

Std  Mean 

0.144494 

Quantiles 

100%  Max 

2.941014 

95% 

2.941014 

75%  Q3 

2.436163 

90% 

2.793092 

50%  Med 

1.96901 

1 10% 

1.230449 

25%  Q1 

1.477158 

5% 

0.819544 

0%  Min 

0.819544 

1 
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Table  9  (Concluded) 


Lewis  and  Clark  sauger  In  seines;  Variable  *  log  (number  per  haul  *  1) 


N 

Mean 

SW  Dev 

Skewness 

CV 

18 

0.596059 

0.35612 

0.888521 

59.74586 

Sum 

Varianoe 

Kurtosis 

Std  Mean 

Quantiles 

100%  Max 

1.491362 

95% 

75%  03 

0.778151 

90% 

50%  Med 

0.477121 

10% 

25%  Q1 

0.30103 

5% 

0%  Min 

0 

Lewie  and  Clark  yeUow  perch  in  seines;  Variable  >  log  (number  pei 

N 

18 

Mean 

0.612634 

Sum 

StdOev 

0.428944 

Variance 

Skewness 

0.488682 

Kurtosis 

CV 

70.01629 

Std  Mean 

Quantiles 

100%  Max 

1.556303 

95% 

75%  Q3 

0.845098 

90% 

50%  Med 

0.60206 

10% 

25%  Q1 

0.30103 

5% 

0%  Min 

0 

10.72906 

0.126822 

1.106837 


1.491362 

1.113943 

0.278754 

0 


11.02741 

0.18399 

-0.27663 

0.101103 


1.556303 

1.146128 

0 

0 


Table  10 

Regression  statistics  by  reservoir,  species,  and  gear 


Fort  Pock  Northern  Pike  tn  Sotnoo;  log  (number/haul  ♦  1) 
Stap  1  Variable  CA6_9  Entered  r-square  -  0.34351 906 


Regression 

Error 

Total 


Variable 

INTERCEPT 

CA6_9 


Estimate 

0.14285371 

0.00004365 


Sum  of  Squares  Mean  Square 
0.33772084  0.33772084 

0.64540026  0.04302668 

0.98312110 


Standard 

Error 

0.05193558 

0.00001558 


Type  II 

Sum  of  Squares  F 
0.32552992  7.57 

0.33772084  7.85 


Prob>F 

0.0134 


Prob>F 

0.0149 

0.0134 


Fort  Peek  Sauger  In  Seines;  log  (number/haul  ♦  1) 

Step  2  Variable  CASUSP  Entered  R-square -0.61910543 


Regression 

Error 

Total 


Variable 
INTERCEPT 
CASUSP 
log(XPA6  94-1) 


Sum  of  Squares  Mean  Square 
0.03138146  0.01569073 

0.01930693  0.00175518 


Parameter 

Estimate 

-13.24362802 

•0.00001168 

3.19689265 


Standard 

Error 

3.17481270 

0.00000448 

0.75848420 


Typed 

Sum  of  Squares 


0.01192039 

0.03118054 


Fort  Peck  Yellow  Perch  In  Seines;  log  (number/haul  ♦  1) 

Step  2  Variable  CASUSU2  Entered  R-square  -0.72756140 


Regression 

Error 

Total 


Variable 

INTERCEPT 

CASUSU2 

log(XPA6_9+1) 


Estimate 

-172.51437506 

-0.00006448 

41.64612989 


Sum  of  Squares 
4.16959679 
1.56132406 
5.73092085 


Standard 

Error 

30.67724333 

0.00002136 

7.33268038 


Mean  Square 

2.08479840 

0.11152315 


Typed 

Sum  of  Squares 
3.52681706 
1.01601045 
3.59740062 


Prob>F 

0.0016 

0.0244 

0.0014 


Prob>F 

0.0001 


Sakakamea  Walleye  In  Gill  Nats;  log  (number/hr) 

Step  1  Variable  CA4_6  Entered  r-square -0.57641 432  C(p)  -  8.37646892 


Regression 

Error 

Total 


Variable 

INTERCEPT 

CA4_6 


Sum  of  Squares  Mean  Square  F 

0.00737146  0.00737146  23.13 

0.00541702  0.00031865 

0.01278848 


Parameter 

Estimate 

0.00760743 


Standard 

Error 


0.00000064 


Typed 

Sum  of  Squares 

0.00050165 

0.00737146 


Prob>F 

0.0002 


Prob>F 

02266 

0.0002 


'independent  variable  abbreviations  are  in  Table  2  and  sample  statistics  are  in  Tables  3-8.  Other  abbreviations  include 
R-square  -  coefficient  of  determination  (multiple  regression);  r-square  -  coefficient  of  determination  (single-variable 
regression);  OF  *  degrees  of  freedom;  F  -  F  statistic;  Prob>F  -  equation  probability. 


TaHs  10  (Continued) 


Sakakawea  White  Crappla  hi  Frame  Nate;  log  (number/hr) 

Step  4  Variable  LINFV4_6  Removed  R -square  -  0.62382038 

Variable  LMG_WY  Entered 


Regression 

Error 

Total 


Variable 

INTERCEPT 

log(MQ_WY+1) 

log(INFV6_9+1) 

CASUSU 


Parameter 

Estimate 

-0.14042899 

-1.54484945 

1.14759144 

0.00000888 


Sum  of  Squares 
0.09000846 
0.05427740 
0.14428586 


Standard 

Error 

0.58533530 

0.64183917 

0.48833657 

0.00000267 


Type  II 

Sum  of  Squares 

0.00020827 

0.02096272 

0.01998315 

0.03990683 


Prab>F 

0.0017 


Prob>F 

0.8136 

0.0294 

0.0329 

0.0047 


Oehe,  NO,  Walleye  In  Qlll  Nets;  log  (number/hr) 

Step  4  Variable  LMG_WY  Removed  R-square  -  0.71 358759 
Variable  LINFV6_9  Entered 


Regression 

Error 

Total 


Variable 

INTERCEPT 

log(INFV6_9+1) 

CA6_9 

CASUSP 


Parameter 

Estimate 

-1.03527621 

0.25772691 

-0.00000648 


Sum  of  Squares 
0.01962360 
0.00787632 
0.02749991 


Standard 

Error 

0.34446804 

0.08401390 

0.00000210 


Mean  Square 

0.00654120 

0.00060587 


Type  II 

Sum  of  Squares 

0.00547*61 

0.00570160 

0.00579732 

0.01427635 


Prob>F 

0.0101 

0.0090 

0.0086 

0.0003 


Oehe,  ND,  White  Bass  In  Qlll  Nsts;  log  (number/hr) 

Step  3  Variable  CA6_9  Entered  R-square -0.64781514 


Regression 

Error 

Total 


Variable 

INTERCEPT 

log(INFV4_6+1) 

CA6_9 

CASUSP 


DF 

Sum  of  Squares 

3 

0.18933719 

13 

0.10293321 

16 

0.29227040 

Parameter 

Standard 

Estimate 

Error 

•1.90059872 

0.97460668 

0.24223918 

-0.00000973 

0.00001085 

0.00000297 

0.06311240 

0.00791794 


Type  II 

Sum  of  Squares 

0.03011163 

0.03405565 

0.02322475 

0.10586342 


Oahe,  NO,  White  Crapple  In  Frame  Nets;  log  (number/hr) 

Step  5  Variable  CASUSP  Removed  R-square  -  0.52495658 
Variable  CASUSU  Entered 


Regression 

Error 

Total 


Variable 

INTERCEPT 

log(MG_WY+1) 

CA6_9 

CASUSU 


Parameter 

Estimate 

-5.96337033 

2.11101896 

-0.00003262 

0.00001323 


Sum  of  Squares 
0.70914581 
0.64171985 
1.35086565 


Standard 

Error 

3.03804471 

1.03558236 

0.00001687 

0.00000639 


Mean  Square 

0.23638194 

0.04936307 


Type  II 

Sum  of  Squares 

0.19019424 

0.20512429 

0.18465847 

0.21162636 


Prob>F 

0.0184 


Prob>F 

0.0714 

0.0624 

0.0752 

0.0589 
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Table  10  (Continued) 


Oahe,  M),  Yellow  Perch  In  Frame  Nats;  log  (number/hr) 

Step  1  Variable  CASUSU  Entered  r-square  -  0.38701 550 


Regression 

Error 

Total 


Variahla 

INTERCEPT 

CASUSU 


Estimate 

0.19600134 

0.00001363 


Sum  of  Squares 
0.42290301 


1.09272885 


Standard 

Error 

0.05149448 

0.00000443 


Mean  Square 

0.42290301 

0.04485506 


Type  II 

Sum  o(  Squares 

0.66021514 

0.42290301 


Proto>F 

0.0077 


Prot&F 

0.0016 

0.0077 


Oahe,  SO,  Northern  Pike  In  Selnee;  log  (number/haul+1) 

Step  6  Variable  LMG_WY  Removed  R -square  -  0.64852268 

Variable  LINFV6  9  Entered 


Regression 

Error 

Total 


Variable 

INTERCEPT 

log(INFV6_9+1) 

CA6  9 

CASUSU2 

CASUSP2 


Sum  of  Squares  |  Mean  Square 
0.01541332 

.00059668 

0.02376682  | 


Parameter 

Estimate 

-0.92134731 

0.22892196 

-0.00000317 

-0.00000281 


Standard 

Error 

023523221 

0.05812956 

0.00000115 

0.00000126 

0.00000139 


Type  II 

Sum  of  Squares 

0.00915363 

0.00925382 

0.00454586 

0.00294881 

0.00477238 


Oahe,  SO,  Walleye  In  Selnee;  log  (numtoer/hauUI) 
Step  2  Variable  LINFV4_6  Entered  R-aquare  < 


023406153 


DF 

Sum  of  Squares 

Mean  Square 

F 

Regression 

2 

0.14506745 

0.07253373 

2.44 

Error 

16 

0.47471595 

0.02966975 

Total 

18 

0.61978340 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

INTERCEPT 

2.35186019 

1.39860245 

0.06389736 

2.83 

log(INFV4  6+1) 

0.63629816 

025569630 

0.09494596 

320 

CASUSU2 

-0.00000446 

0.00000245 

0.09873120 

323 

Oahe,  SO,  White  E 

lass  In  Seines;  log  (number/haul+1) 

Step  4  Variable  CA6  9  Removed 

R -square  -  0.69067754 

Variable  CASUSU  Entered 

DF 

Sum  of  Squares 

Mean  Square 

F 

Regression 

3 

224654767 

0.74884922 

11.18 

Error 

15 

1.00612458 

0.06707497 

Total 

18 

325267225 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

INTERCEPT 

127989249 

0.07129923 

21.61412998 

32224 

CA4  6 

0.00002130 

0.00000927 

025466639 

529 

CASUSU 

l.l.Y.Y<'<Y:«*fl 

0.00000874 

0.96073432 

14.32 

CASUSU2 

-0.00001831 

0.00000500 

0.89756816 

13.38 

Prob>F 

0.0037 


Prob>F 

0.0015 

0.0015 

0.0153 


0.0134 


ProboF 

0.1185 


Prob>F 

0.1121 

0.0926 


Prob>F 

0.0004 


Prob>F 

0.0001 


0.0018 
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Table  10  (Continued) 


Oahe,  SO,  Whit*  Crapplo  In  Seinee;  log  (number/haul+1) 
Stop  2  Variable  CA6_9  Entered  R -square  -  0.661 1 31 60 


Regression 

Error 

Total 


Sum  of  Squares 
1.182S9234 
0.60614681 
1.78873915 


Mean  Square 

0.59129617 

0.10102447 


ProtoF 


Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

■0.30236343 

0.21707221 

0.19600928 

1.94 

02131 

CA6  9 

•0.00008683 

0.00002873 

9.14 

CASUSP 

0.00005184 

0.00001603 

1.05654080 

10.46 

0.0178 

1  Oarie,  SO,  Yellow  Perch  In  Seines;  log  (number/haul+1) 

|  Step  1  Variable  LMG_WY  Entered 

i- square  -  0.42090188 

OF 

Sum  of  Squares 

Mean  Square 

F 

Prob>F 

Regression 

1 

3.96568961 

3.96568961 

12.36 

0.0027 

Error 

17 

5.45619666 

0.32095274 

Total 

18 

9.42188627 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

•13.30048861 

4.18661985 

3.23928551 

10.09 

log(MG  WY+1) 

5.12792536 

1.45882394 

3.96568961 

12.36 

0.0027 

Sharpe  Qizzard  Shad  In  Seines;  log  (number/haul+1) 

Step  1  Variable  LFR4_6  Entered 

r-square  -  0.46707387 

DF 

Sum  of  Squares 

Mean  Square 

F 

Prob>F 

Regression 

1 

1.67717589 

1.67717589 

14.02 

0.0018 

Error 

16 

1.91363818 

0.11960245 

Total 

17 

3.58081507 

Parameter 

Standard 

Typed 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

4.77853781 

0.66118699 

6.24714369 

0.0001 

log(FR4_6+1) 

-3.94488441 

1.05345268 

1.67717589 

14.02 

0.0018 

Sharpe  Walleye  In  Seines;  log  (number/haul+1) 

Step  1  Variable  LMG_WY  Entered  r-scuare  -  0.4770221 0 


Regression 

Error 

Total 


Variable 
INTERCEPT 
log(MG_WY+1 ) 


Parameter 

Estimate 

17.88154868 

•6.11484906 


Sum  of  Squares 
0.58869534 
0.64540962 
123410495 

Mean  Square 

0.58869534 

0.09220137 

F 

6.38 

Standard 

Typed 

Error 

Sum  of  Squares 

F 

6.78859063 

0.63971776 

6.94 

2.41996531 

0.58869534 

6.38 

Prob>F 


Prob>F 

0.0337 
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Table  10  (Continued) 


8akakawea  Whit*  Crappie  In  Frame  Nats;  log  (numtoer/hr) 

Stop  4  Variable  LINFV4_6  Removed  R-square  -  0.62382038 
Variable  LMV  WY  Entered 


Regression 

Error 

Total 


Variable 

INTERCEPT 

log(MG_WY+1) 

lofl<INFV6_9+1) 

CASUSU 


Parameter 

Estimate 

•0.14042899 

-1.54484945 

1.14759144 

0.00000888 


Sum  of  Squares 
0.09000846 
0.05427740 
0.14428586 


Standard 

Error 


0.64183917 


0.00000267 


Mean  Square 


0.00361849 


Type  II 

Sum  of  Squares 
0.00020827 
0.02096272 
0.01998315 


Prot»F 

0.0017 


Prob>F 

0.8136 

0.0294 


0.0047 


Oatie,  NO,  Walleye  In  Qlll  Nets;  log  (number/hr) 

Step 4  Variable LMG.WY Removed  R-square  - 0.71 358759 
Variable  LINFV6  9  Entered 


Regression 

Error 

Total 


VarWrle 

INTERCEPT 

bfl(INFV6__9+1) 

CA69 

CASUSP 


Oalte,  NO,  White  Bass  In  Gill  Nets;  log  (numbor/hr) 

Step  3  Variable  CA6_9  Entered  R-square  >0.64781514 


OF 

Sum  of  Squares 

Mean  Square 

3 

0.01962360 

0.00654120 

13 

16 

0.00787632 

0.02749991 

0.00060587 

Parameter 

Standard 

Type  II 

Estimate 

Error 

Sum  of  Squares 

-1.03527621 

0.34446804 

0.00547261 

0.25772691 

0.08401390 

0.00570160 

-0.00000648 

0.00000210 

0.00579732 

0.00000389 

0.00000080 

0.01427635 

Oahe,  ND,  White  Crapple  In  Prams  Nets;  log  (number/hr) 

Step  5  Variable  CASUSP  Removed  R-square  -  0.52495658 

Variable  CASUSU  Entered 


Prob>F 

0.0008 


Prob>F 

0.0101 

0.0090 


Regression 

Error 

Total 

DF 

3 

13 

16 

Sum  of  Squares 
0.18933719 
0.10293321 
029227040 

Mean  Square 

0.06311240 

0.00791794 

F 

7.97 

Prot»F 

0.0029 

Variable 

Parameter 

Estimate 

Standard 

Error 

Type  II 

Sum  of  Squares 

F 

Prot»F 

INTERCEPT 

-1 .90059872 

0.97460668 

0.03011163 

3.80 

0.0731 

k>g(INFV4  6+1) 

0.50238066 

0.24223918 

0.03405565 

4.30 

0.0585 

CAS  9 

•0.00000973 

0.02322475 

2.92 

0.1105 

CASUSP 

0.00001085 

0.00000297 

0.10586342 

13.37 

0.0029 

Regression 

Error 

Total 

DF 

3 

13 

16 

Sum  of  Squares 
0.70914581 
0.64171985 
1.35086565 

Mean  Square 

023638194 

0.04936307 

F 

4.79 

Prob>F 

0.0184 

Variable 

Parameter 

Estimate 

Standard 

Error 

Typo  II 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

-5.96337033 

3.03804471 

0.19019424 

3.85 

0.0714 

log(MG  WY+1) 

2.11101896 

1.03558236 

020512429 

4.16 

0.0624 

CA6  9 

0.00003262 

0.00001687 

0.16465847 

3.74 

0.0752 

CASUSU 

0.00001323 

0.00000639 

021162636 

429 

0.0589 
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Table  10  (Continued) 


Sharp*  Whit*  Bass  In  Seines;  log  (number/haul+1) 

|  Step  1  Variable  LFR-i  6  Entered 

r-square  -  0.33342278 

DF 

Sum  of  Squares 

Mean  Square 

F 

Prob»F 

Regression 

1 

1.42944552 

1.42944552 

8.00 

0.0121 

Error 

16 

2.85774063 

0.17860879 

Total 

17 

4.28718615 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

2.90216571 

0.80798942 

230431516 

12.90 

0.0024 

iog(FR4  6+1) 

-3.64190641 

1 28734927 

1.42944552 

8.00 

0.0121 

1  Sharp*  Freshwater  Drum  In  Seine*;  log  (numfaer/hauUI 

) 

|  Step  1  Variable  LINFV4_6  Entered 

r-square  -  0.53057930 

DF 

Sum  of  Squares 

Mean  Square 

F 

Prob>F 

Regression 

1 

1.01830204 

1.01830204 

18.08 

0.0006 

Error 

16 

0.90092480 

0.05630780 

Total 

17 

1.91922683 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

9.44482884 

2.12292895 

1.11451475 

19.79 

0.0004 

log(INFV4  6+1) 

-2.35464544 

0.55369614 

1.01830204 

_ 

18.08 

0.0006 

1  Francis  Case  Gizzard  Shad  in  Seines;  log  (number/haul+1) 

|  Step  1  Variable  CA6_9  Entered 

r-square  -  0.41412201 

DF 

Sum  of  Squares 

Mean  Square 

F 

ProboF 

Regression 

1 

1.11271293 

1.11271293 

4.95 

0.0615 

Error 

7 

1.57420759 

022488680 

Total 

8 

2.68692052 

Parameter 

Standard 

Typed 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prot»F 

INTERCEPT 

2.56307236 

0.31066960 

15.30693327 

68.07 

0.0001 

CA6_9 

0.00032631 

0.00014670 

1.11271293 

4.95 

0.0615 

1  Francis  Case  Walleye  In  Gill  Nets;  log  (number/net  night 

) 

|  Step  2  Variable  LINFV4_6  Entered 

R-square  -0.939081 13 

DF 

Sum  of  Squares 

Mean  Square 

F 

Frob>F 

Regression 

2 

0.61435592 

0.3071796 

4625 

0.0002 

Error 

6 

0.03985371 

0.00664229 

Total 

8 

0.65420963 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Proo>F 

INTERCEPT 

-12.75772262 

2.08084273 

024968092 

37.59 

0.0009 

log(INFV4  6+1) 

3.22437950 

0.53245589 

024358054 

36.67 

0.0009 

log(SBINF46+1) 

0.66742812 

0.07126818 

0.58255271 

87.70 

0.0001 
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Table  10  (Continued) 


Francis  Csss  White  Bass  In  Sainss;  log  (number/hauM) 
Stop  1  Variable  CA4_5  Entered  r-equare  -  0.46677001 


— 

Regression 

Error 

Total 

DF 

1 

7 

8 

Sum  of  Squares 
0.89064810 
1.01746101 
1.90810911 

Mean  Square 

0.89064810 

0.14535157 

F 

6.13 

Prob>F 

0.0425 

Variable 

Parameter 

Estimate 

Standard 

Error 

Type  II 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

1.50305607 

0.18437538 

9.65971309 

66.46 

0.0001 

CM  5 

0.00026772 

0.00010815 

0.89064810 

6.13 

0.0425 

Francis  Casa  White  Crappte  In  Sainss;  log  (number/haul+1) 
Step  1  Variable  CASUSP2  Entered  r-aquare  -  0.3650461 1 


DF  Sum  of  Squares  Mean  Square  F  Prob>F 

Regression  1  1.90714948  1.90714949  12.41  0.0023 

Error  19  2.92051352  0.15371124 

Total  20  4.82766301 


Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

0.14198927 

0.09977605 

0.31128918 

2.03 

0.1709 

CASUSP2 

0.00012352 

0.00003507 

1.90714949 

12.41 

|  Francis  Case  Yellow  Porch  In  Seines;  log  (number/hauM) 

Step  1  Variable  CM_5  Entered 

r -square  -  0.85184871 

DF 

Sum  of  Squares 

Mean  Square 

F 

ProtfeF 

Regression 

1 

0.94177748 

0.94177749 

40.28 

0.0004 

Error 

7 

0.16379147 

0.02339878 

Total 

8 

1.10556896 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

0.78560036 

0.07397580 

2.63886524 

112.78 

0.0001 

CM  5 

0.00027530 

0.94177749 

4025 

0.0004 

Lewis  and  Clark  Gizzard  Shad  in  Seines;  log  (number/hauM) 

Step  2  Variable  LFR4_6  Entered 

R-square  -  0.28563806  C(p) -4.79250172 

DF 

Sum  of  Squares 

Mean  Square 

F 

Prot»F 

Regression 

2 

1.82488060 

0.91244030 

3.00 

0.0802 

Error 

15 

4.56390599 

0.30426040 

Total 

17 

6.38878659 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

9.55360822 

3.41755804 

2.37765309 

7.81 

0.0136 

LSBINF46 

-2.23951226 

0.96401494 

1.64204400 

5.40 

0.0346 

LFR4  6 

-2.55878867 

1.93304204 

0.53312822 

1.75 

02054 

(Shmt7of8)\ 


Table  10  (Concluded) 


Lewis  and  dark  Yellow  Parch  in  Seines;  log  (number/haul+1) 

|  Step  2  Variable  LINFV4_6  Entered 

R-square  -  0.32665666 

DF 

Sum  of  Squares 

Mean  Square 

F 

Prob>F 

Regression 

2 

1.02174163 

0.51087081 

3.64 

0.0515 

Error 

15 

2.10613467 

0.14040898 

Total 

17 

3.12787630 

Parameter 

Standard 

Type  II 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prot»F 

INTERCEPT 

8.00092271 

5.69663642 

0.27697338 

1.97 

0.1805 

iog(INFV4  6+1) 

-0.99428531 

1.32519435 

0.07904197 

0.56 

0.4647 

iog(SBINF46+i) 

-1.79395113 

0.66518762 

1.02123992 

7.26 

0.0166 

|  Lewie  and  Ctark  Sauger  In  Seines;  log  (nurnber/haul+1) 

1  Step  2  Variable  LFR6_9  Entered 

R-square  -  0.49603852 

Regression 

DF 

Sum  of  Squares 

Mean  Square 

F 

Prub>F 

2 

1.06944426 

053472213 

7.38 

0.0059 

Error 

15 

1.06652591 

0.07243506 

Total 

17 

2.15597017 

Parameter 

Standard 

Typed 

Variable 

Estimate 

Error 

Sum  of  Squares 

F 

Prob>F 

INTERCEPT 

7.19505456 

1.86541787 

1.07761778 

14.88 

0.0016 

LSBINF46 

-1.62353292 

0.46797126 

0.87183173 

12.04 

0.0034 

LFR6  9 

-2.45795286 

1.02712074 

0.41481355 

5.73 

0.0302 

(Sh—t8otB) 


Table  11 

Percent  of  years,  under  the  large-seasonal-drawdown  alternative.  In  which  the 
value  of  an  independent  variable  was  1, 10, 20, 30,  SO,  or  100  percent  outside  the 
range  of  the  original  data  used  to  derive  regression  equations 


FrandsCaae 

CASUSP2 

SBINF 

CA4  5 

INFV4  6 
CA6_9 

Fort  Pack 

CASUSP 
CASUSU2 
CA6  9 
XPA6_9 

Lewis  and 
Clark 

SBINF 

INFV4  6 

FR4  6 

FR6_9 

Oahe 

MG  WY 

CASUSP 

CASUSU 

CASUSP2 

CASUSU2 

CA4  6 

INFV4  6 

CA6  9 
INFV6_9 

Sakakawea 

MG  WY 
CASUSU 

CA4  6 
INFV6_9 

Sharpe 

MG  WY 

INFV4  6 

FR4  6 
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